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Abstract 

The Dirac method of quantizing Hamiltonian systems with constraints is 
applied to the massless Thirring model. We solve the quantum Hamiltonian 
equation for the energy-momentum tensor and obtain a violation of the classical 
conservation law. A previously noticed problem with the equal-time anticommu- 
tators can be fixed using this Hamiltonian method. 
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1. Introduction 

Since the late 1960's an extensive literature has evolved on the massless 
Thirring model. We would like to fill a gap, however. 

The Thirring model is a typical system with constraints and so we quantize 
it here according to Dirac's special Hamiltonian formalism. Unlike most authors, 
however, we don't use solutions of the normal-ordered Lagrange equations [1-4] in 
our treatment, since such solutions are not part of the Dirac-Hamiltonian formal- 
ism [5]. We only consider operators, commutators, and normal ordering in initial 
time-like elements. Only Hamiltonian language is used and the quantum Hamil- 
tonian equations of motion are solved. We don't mix elements of the Lagrange 
and Hamiltonian formalisms. 

This method was successfully tested in the quantum integrable sine-Gordon, 
Zhiber-Shabat, nonlinear Schrodinger, Korteweg-de Vries and modified Korteweg- 
de Vries systems [6,7]. Quantum commutative integrals of motion for these mod- 
els were constructed as solutions of the quantum Hamiltonian equations. This 
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is the unique known nonperturbative way of constructing quantum integrals of 
motion, at the present time [6]. In addition, using this method in the massless 
Thirring model helps to reveal a new infinite dimensional symmetry in the sine- 
Gordon theory, explaining the power law behavior of its correlation functions 
[7]. 

Furthermore using this Hamiltonian approach in the massless Thirring model 
we avoid a problem with equal-time anticommutators. In fact we exchange this 
problem in its hidden form (in the equal-time anticommutator) for an obvious 
anomaly. Indeed the difference between the standard equal-time anticommutator 
from the canonical one is of order g'^ for {g <^ 1). We fix this problem in the 
anticommutator, but get instead a correction in the same order of g in another 
place: in the degree of non-analyticity of the operator That is the main 

result of this article. It is important because this anomaly breaks the well known 
conformal symmetry of the massless Thirring model. 

In Sec. 2 we consider the naive solution of this model in a modern, adapted 
version [3] and explain the unsolved problem with this approach. In Sec. 3 we 
use the classical Hamiltonian formalism for systems with constraints on the mass- 
less Thirring model. Wc solve the simple constraints and calculate their Dirac 
brackets. More complicated constraints have complicated Dirac brackets and so 
we work with them in a special way: very useful operators P± are introduced. 
In Sec. 4 we establish an important property of initial functions, and so we use 
these functions in all operators of our theory. In Sec. 5 we calculate a quantum 
anomaly of the components of the energy-momentum tensor, a quantum correc- 
tion to the classical conservation law. Sec. 6 is devoted to concluding remarks and 
a brief discussion of the connection of this work with the Thirring/sine- Gordon 
equivalence. 



2. Problem with equal-time anticommutator 

Let us recall aspects of the massless Thirring model in the modern conformal 
field theory approach. We have the vertex operator [3] 

Vm,n{z,z) =: exp2z r+ ^ ^ e{z) + ^+ ^ ^ £(^)j : (2.1) 

so that e{z,z) satisfies equation d^dze^z^z) = 0, where e{z,z) = e{z) + e{z). 
Components of the energy- momentum tensor are given by T{z) ~: (e(z))^ : and 
T{z) ~: {e{z)y :. In the standard approach these are holomorphic and antiholo- 
morphic quantities. Conformal dimensions have the form 

{A,A) = {P+m + p.n)M(3+m-P-nf), P+P- ^ ^, (2.2) 
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We can consider the obvious equations 

dzVm,niz, z) =: i{(3+m - (3-n)d^e{z))Vm,n ■, 

{z, z) =: i{(3+m + (3-n)dze{z))Vm,n '■ ■ (2.4) 
After identifying certain vertex operators with fields of theory 

= iJ2 = V._^ M-V.., V'^^V^ii, (2.5) 

we have equations for massless Thirring model 

dzi^i = g ■ Ji>i dziJ2 = g ■ J^i ■ ■ (2.6) 
These solutions for ■0i,V'2 satisfy normal ordered Lagrange equations. 

But for a quantum solution we also must demand the correct equal-time 
anticommutator 

[Mx, t),i;t{y, t)]+ = iS{x - y). (2.7) 

This anticommutator is a fundamental property of the quantum theory. In the 
considered case above the operator algebra of the solutions has the form 

'4)t{z,z)i)i{Q,Q) = z ^ . (2.8) 

In fact it is possible to prove the cqual-timc property (which follows from the 
above) only for the (/ = case. And so this solution has a problem, as indeed we 
will show the quantum solution does. 

Let us consider one example of an explanation in the literature [4] . The author 
considers the solution 

'(/'(x) = exp[— z67^0*-^-'] exp[m0^~''(x)] exp[m0*-^-*(a;)] e'iq)[—ilr)^4'^^\x)]u (2.9) 

of the Thirring model, where m is a "two component c-number quantity". For 
the equal-time anticommutator 

< 0|[V'.(x), V';(?/)]+|0 >= ur{x)u^{y)5{x - y)(x - yf-^ (2.10) 

is obtained. Here h — ^-^^ > 1, « and h are constants in the theory. By requiring 

Ur{x)uUy) {x - y)-^^-'\ (2.11) 

the author gets the ordinary 5 function result for the cqual-timc anticommutator. 
But we believe that it is not possible to construct this "c-number quantity" in a 
space of functions. Indeed we have from (2.11) 

\ogUr{x) + log^^+(y) ^{l-h) \og{x - y), (2.12) 
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(for noncommutative operators a,b; we have log(a6) ^ logo + log 6). If dxdy 
operator acts on the left side, (2.12) we get 0; on the right side we get -^E^ 0- 

In the Hamiltonian approach, the canonical equal-time relation is postulated 
at the beginning of the calculation and is preserved in time evolution, and so this 
problem disappears. 

If we consider the massive Thirring model we get a similar problem with the 
equal-time relation for fermions when 5^ 7^ [8] 

[i;{x, t),ij+iy, t)]+ = ix- yr^^'^6{x - y) ^ 5{x - y). (2.13) 

This problem can also be fixed by the Hamiltonian method of quantization. We 
will consider the massive case in a separate work. 



3. Classical Hamiltonian formalism 

We use light cone coordinates with the notation 

x± = i=(x°±xi), (3.1) 
and represent the 7^-matrices by 

If we think of x^ as the imaginary time ir, we work in Euclidean space, and 
x"^ ~ ^, ~ z are its complex coordinates. 

We start with the Lagrangian for the massless Thirring model in light cone 
coordinates 

>C = ^ {^td+i^2 + - - - 2gi,t^i'4^t'4^2 . (3.3) 

The canonical conjugates of the fields are the following: 

^'^^ ^ ^ "^"^^ ' ^ '''''' 

dL iV2 , ,1 iV2 , , 

- _ n f 1 - 
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We must have the canonical Poisson brackets (more exactly Poisson-Berezin 
brackets for our anti-commutative variables) for the fields and their conjugate 
fields: 

{i(ji{x);7r^,{y)} = 5{x-y), {'ilJi{x);7T^^{y)} = 5{x-y), 
{ilJ2{x);7r^2{y)} = 6{x -y), {^^ (a;); 7r^+ (2/)} = - y). (3.5) 

The expressions = 0, (i = 1, . . . , 4) are the primary constraints. We have the 
first step Hamiltonian density, 

A2 (-^^1 + T^v^) + ^^37^^-2 + A47r^+ , (3.6) 

constructed in the usual way for systems with constraints [5] . We must demand 
the conservation of constraints in time (a dot indicates a derivative with respect 
to "time" x_). From fl = ff = 0, wc can obtain A];,A^, and from — fl — 0, 
we find new (secondary) constraints: /|, /| 

fl = {H,- fl} = 0; A} = -^^1^2+ 

fl = {i/i; fl} = 0; A^ = _-^V'i+V'2+V'2, 

fl = {H,- fl} = 0; fl = tV2d+ij2 - 2g^ti^,^2 = 0, (3.7) 

fl = {Hi; fl} = 0; fl = zV2d+i;+ + 2#+ViV'2^ = 0. 

Here we have introduced Hi = J hidz. The constraints fl,fi are part of the 
Lagrange equations but in the Hamiltonian sense, they are only constraints of 
second class. We must introduce the second step Hamiltonian density 7^2: 

-2gijti^t^2 (^^1 + v) + ^37^^2 + Al7r^++ (3.8) 

+A^ (iV2d+ii^2 - 2t/^+Vi'02) + xl (iV2d+i^t + '^gi^ti^ii^t) ■ 

Demanding the conservation of the new constraints in time yields 

fl^{H2;fl}^0; Xl 
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f!^{H2-j!}^0; Xl, (3.9) 

fl = {H,-fl} = 0; \l = 0, 

fl = {H2;fl} = 0; \l = 0. 

From /I = /I = we can determine A3, A4, but their forms are not important for 
this work. Similarly, /a = /i = determine A|, A4. Thus all the constants in the 
H2 Hamiltonian density (3.8) are determined and we have no new constraints. 

It is very useful to resolve the constraints fl = = ^ and so we must 
calculate Dirac brackets. Using the notation (/i = fi;f2 — /2,Q;,/? = 1,2) the 
expressions for the Dirac brackets are [5] 

{i^l, Ipt joirac = {i^l,1pt}- I^iV'l, fa}{fa, Uy^{fa, Ipt }, 
{^l(a;),^/;i(y)}Dirac = {V'l(a;)"^,^i^(l/)}Dirac = 0. 

In our case, the matrix of the constraints {fafp} is not degenerate, and so has 
an inverse. We used 5~^{x — y) = 5{x — y). In our theory the physical anti- 
commutative variables are ipi, ipf only, and the anti-commutative variables -021 i^t 
are dependent. We don't resolve the f^, fl, /|, /| constraints, and so we will 
use the ordinary Poisson (3.5) brackets between the fields 1^2, 4^2+ and Tr^a^'^vJ' 
but then impose the constraints. 

After resolving the constraints fl = fl = 0, the Hamiltonian density 7^2 has 
the form: 

In our research, we are only interested in functionals I[i(ji,'iIjI~]. The important 
part of 7-^2 is therefore 7^2 = '^g'4'ii'i'4'2'4'2, since {7^2, = {'^21 ^bPi^i'i]}- 

An important remark can now be made. The constraints /|,/| must be 
"hamiltonized" . We must introduce fields P± as integrals over certain densities: 

—^J i^t^i'^27r^2dz , P-^ -—^J i/jfi^ii^t^i^+dz . (3.12) 

Here tt^j and 7r^+ are our constraints fl and fl- The operators P± do not vanish, 
however, because their action is defined so that the constraints are imposed only 
after calculating the Poisson brackets. 
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After quantization, the operators P± help to remove singularities in the theory. 
We calculate 

d+1p2 = {P+, -02} = J Ipt iJl1p2T^i,2dz, ^2} = J^^l (3.13) 

using {7/^2(3;), 7r^2(y)} = ^{.^~y)-i ^'^'^ obtain the constraint /|. Another important 
example of the action of these operators is 

2g f 

d+{ll^2)A = {P+,1p2iJl} = — /|{ J i^t A'ip2T^^2dz,1p2'ipi} = 

Vl" V'lV^2V'l + #lV'27rv,2V'2 = —j^'iPi i'ii^2i>u (3.14) 



,,,,,, , 
where we used 

{V'l(x),'0]^(y)}Dirac ^^5{x-y), {V^l(a;), '0l(?/)}Dirac = 0, 

{'il)2{x),n^^{y)} = 5{x - y), = 0. (3.15) 

While d+ acts on -02 only, P+ acts on all fields 7/^2 and ■(/'i, and imposing the 
constraints removes the extra parts in the classical case. We will consider the 
quantum analog of this example, and will find a nontrivial action of this operator. 
Expression (3.14) is equal to zero because ipj — 0. We must remember here that 
constraints n^^ — ~ must be imposed after calculating Poisson brackets 
(or commutators (anti-commutators) in the quantum case). 



4. Quantization 

Recall the solution of the quantum Hamilton equation 

il^{z,z) = exp (^—izH^ 'i^'^{z,z) exp (izH^ . (4.1) 

The superscript indicates an initial quantum field, and ijj{z, z) denotes an op- 
erator solution. Let us consider the operator Hamilton equation of motion: 

i^J{z,z)^[H,^P{z,z)] (4.2) 

Inserting the solution gives 

id^i^iz, z) = [H, tpiz, z)] + i exp (^-izH) d^ip^iz, z) exp (izH) (4.3) 
For consistency then, we must demand that the initial operator obey 

d-'ii^\z,z) = 0. (4.4) 
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This is a very useful property (that of analytical functions) for quantization. 

It is a simple exercise to check that the solution above satisfies fermionic prop- 
erties too. For example, il^^{z,z) — (this follows from the fermionic properties 
of the initial quantum fields). It is easy to check property in Hamiltonian ap- 
proach [iIj{z, z), iI^~^{z', z)]+ = id{z — z') (we have postulated for initial operators 
[il^^{z,z),il^^~^{z',z)]+ = iS{z — z')). And so we don't have the problems like in 
naive approach. 

We will use the initial fields in all operators (H, P±) and will drop the "0" and 
z in the notation '^^(z,^). The Poisson brackets above (for the classical initial 
fields) can be quantized in the usual way, and we obtain the standard singular 
parts of operator product expansions: 



{z - z') ' 



1 

{z~z') ' v-y^jv-/ - (^_^/)- 



Here we used the standard expansions of the analytical fields: 

n n 

^2(0 = E^2(e-«)-"-^ ^2+(0 = EV'2+'^(e-«)-"-^ (4.6), 

n n 
n n 

where a is the center of the expansion. We have also introduced some redetermine 
all of the fields multiply by unimportant constants. 

Poisson brackets for the complicated operators H, P± must also be quantized. 
The integration contours for both parts of the commutator (Q — H, P±): 

[Q, A{i')] = QA{i') - A{i')Q, Q = J q{Od^- (4-7) 

need to be determined. Recall also that we are interested in only a certain part 
of the Hamiltonian, Our quantum operators are 

P+ = — ^ J ■ '(Pt V'iV'27rv,2 -.dx, P_ = -— ^ J : tpt ipiij^ Tt^+ : , (4.8) 
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where :: denotes normal ordering at the initial time. Let us consider the contour 
integration in our commutators. In the first part of the commutator, we choose 
the contour closing above the point ^' (see Fig.l). 

In the second part, we choose the contour closing below the same point ^' (see 
Fig.2). 

Fig 2. 




Fig. 2 Contour of the integration for the commutator. 



Letting the radius R of the semicircles go to oo, we have 

/ n',{OA{ad^^0, f p±{i)A{Odi^Q, (4.9) 

using the asymptotic behavior of the fundamental fields ~ So for our 

operators H2, P±, (using analytic property (4.4)): 

[Q,MO] = f^jmim, (4.10) 

where the notation indicates closed-contour integration around the point ^ = 

Let us consider the action of the quantum operators P± (4.8). In our case we 
must calculate expressions like d+{ip2)'ipi- The naive way (without introducing 
P± operators) is would be 

lim,a+(t^2)(^)^i+(;^') = ^ lim/^+^i^2)(^)^i+(^')- (4.11) 

We find a singularity for z — > z', arising from the product il}i{z)ipi {z') ~ ^3^, so 
that 

lim{i>ti'ii'2){z)i'tiz') = oo- (4-12) 
9 



If we use the action of the quantum operators P±, we instead obtain 

d+{iP2)i^t = whd+{iPt'4^2). (4.13) 

The singularity has been removed, but there is a quantum correction to the 
classical result [h is Planck's constant). This is the idea behind the introduction of 
these operators P±. Using the quantum Hamiltonian H, when we calculate [H, I] 
we find singularities. The analytic property of initial fields (4.4), however, with 
the choice of the contour of the integration described above help to remove those 
singularities. Incidentally, this way (in [H,I]) of removing the singularity gives 
the correct commutative integrals of motion (elements of Hamiltonian formalism 
too) for quantum sine-Gordon theory [7]. 

In this article we will solve the quantum Hamilton equation 

id_A - [H,A], (4.14) 

for the "++" component of the energy-momentum tensor. That is, we put A = 
T+_|_ = T, with 

f=:i^id+ijt ■ + -i^td+i;i:, (4.15) 
where tpiji^t ^-^^ T are initial functions. We must calculate [H,T]. 



5. Calculation of anomaly 

Let us consider 

i2 =: ijjid+ipf : . (5.1) 
After a simple calculation we obtain 

[hM- (5.2) 

= 2gi{-) : i,Jti'2dJi : +2gt{-) : 9+ (^'+4+ ^■2)^1+ : +gthdl{^t^2){-). 
Similarly, for 

4 = - : V'^^+V'i (5.3) 

we find 

(5.4) 

- 2gi{-) : ^1^+^25+^1 : +2gt{-) : d+{^t^t^2)^t ■ +gthdl{tP+iJ2){-). 
The result is 

[H,f] = Agi:iP,d+iijt^2)ijt-- (5.5) 
Now we can calculate : dj^{il)2^2 ) '■ using our operators P±. We find 

(5.6) 
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+ {-igh')dl^ti'2 , 

+ {-igh')dli,2'^t- 

We therefore have 

: -V'i^i+a+(V'24+) : = -igh{dli^2il^+ + a^V'2+V'2) (5.7) 

yielding 

[H,f] = -2fx\dliP2i't + 5+4+4) , t^-igh. (5.8) 

If [H,T] is not zero, then we have a quantum anomaly, d^T ^ 0. This is 
easily established: 

(1 - ^)(a^44+ + 5+4'"4) = 2z^(: a+4+444+ : + : 4+9+444'' 7^ 0. 

(5.9) 

The "++" component of the energy-momentum tensor (4.15) is therefore not 
conserved in time. However, if we have h — ^ (classical limit) or = (free 
massless fermions), we do have 

[^,f]=0, (5.10) 

and so T is conserve for a quantum free massless fermion theory. 

In the usual quantum case only one mode (momentum) is conserved, and 
indeed we can make a simple transformation to obtain 

\H,f\ = 2A|2a+(a+V'24+ + 5+4+4)- (5.11) 

If we introduce notation for the momentum operator I2 — J Tdz, we have [H, 12] — 

0. If we want to calculate T component of the energy-momentum tensor we 

must consider variable like the time, from the beginning of the calculation. 

6. Conclusion 

An important conclusion can now be drawn. We can bosonize the initial 
functions using 

4(2) = : exp (0(2)) : , '0+(z) : exp (-4>{z)) :, (6.1) 

where is the initial bosonic function (operator) . After a simple transformation 
we find 

f = : 45+4+ : + : 4+^+4 ■ ^ - ■ {d+4>f ■ ■ (6.2) 
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We see the equivalence between the "++" component of the energy-momentum 
tensor for bosonic and Thirring fermionic theories at the initial time. A similar 
calculation can be found in [9]. In the quantum massless Thirring model we 
have f^-Trhirring 7^ 0, however. So, notwithstanding the equivalence of the initial 
operators T, the equivalence between the massless Thirring model and the mass- 
less free bosonic field is lost, if the Hamiltonian formalism is used throughout. 
Indeed let us consider the free massless bosonic theory. It has the Hamiltonian 
-^boTon = in hght cone coordinates, and so T^oson (6-2) is conserved in time, 

Q rpfree _ n 
^ — ^ boson 

We see then, the violation of the zeroth-order approximation (in the sense 
of [10]) of the sine-Gordon /Thirring model equivalence [8,10]. Of course, the 
solution [1-4] can still sometimes be useful, because the violation is very weak, 
~ g^h^, for 5f < 1. 
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